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We report the first observation of a charm-strange meson D^(2632) at a mass of 2632.5 ± 
1.7 MeV/c'^ in data from SELEX, the charm hadro-production experiment E781 at Fermilab. This 
state is seen in two decay modes, D^77 and D^'K"*". In the D^t; decay mode we observe a peak with 
101 events over a combinatoric background of 54.9 events at a mass of 2635.4 ± 3.3 MeV/c^ . There 
is a corresponding peak of 21 events over a background of 6.9 at 2631.5 ± 2.0 MeV/c^ in the decay 
mode D^K"*". The decay width of this state is < 17 MeV/c^ at 90% confidence level. The relative 
branching ratio r(D°K"'")/r(D^r7 ) is 0.14 ± 0.06. The mechanism which keeps this state narrow is 
unclear. Its decay pattern is also unusual, being dominated by the D^t; decay mode. 

PACS numbers: 14.40. Lb, 13.30.Eg, 13.25Ft 

In 2003 the BaBar collaboration reported the first Dj}(2317). A number of theory papers suggested dif- 

observation of a massive, narrow charm-strange meson ferent explanations for the unexpectedly low mass of 

Dj^(2317) below the DK threshold [ij. Confirmation the state, which had been thought to lie above the DK 

quickly followed from CLEO and BELLE 1\. The threshold 0, H H B • A prediction relevant to this 

CLEO collaboration showed that a higher-lying state, experiment was that the pattern of parity-doubled states 

suggested by BaBar, existed and was a partner to the is expected to continue to higher excitations with similar 



splittings |l(|. 

The SELEX experiment used the Fermilab charged 
hyperon beam at 600GeV/c to produce charmed par- 
ticles in a set of thin foil targets of Cu or diamond. 
The negative beam composition was approximately half 
E~ and half 7r~. The three-stage magnetic spectrom- 
eter is shown elsewhere [llj, [l^- The most important 
features are the high-precision, highly redundant, ver- 
tex detector that provides an average proper time reso- 
lution of 20 fs for charm decays, a 10 m long Ring-Imaging 
Cerenkov (RICH) detector that separates tt from K up 
to 165GeV/c 13], and a high-resolution tracking sys- 
tem that has momentum resolution of ap/p< 1% for a 
150GeV/c track. Photons are detected in 3 lead glass 
photon detectors, one following each spectrometer mag- 
net. The photon angular coverage in the center of mass 
typically exceeds 27r. For this analysis the photon en- 
ergy threshold was 2 GeV. Previous SELEX Dg stud- 
ies showed that most of the signal came from the S~ 
beam 'l4!|. We restrict ourselves in this analysis to the 
10 X 10^ E^-induced interactions. 

In this study we began with the SELEX 
Dg — > K''"K~7r sarnple used in lifetime and hadropro- 
duction studies [lj,lla|- Charged conjugate final states 
are included here and throughout this paper. The 
sample-defining cuts are defined in the references. The 
Dg meson momentum vector had to point back to 
the primary vertex with x^ < 8 and its decay point 
must have a vertex separation significance of at least 
8 from the primary. Tracks which traversed the RICH 
(p > 22 GeV/c) were identified as kaons if this hypothesis 
was most likely. The pion was required to be RICH- 
identified if it went into its acceptance. There are 
544 ± 29 S^ induced signal events with these cuts. 

Due to high multiplicity, photon detection in an 
open charm trigger is challenging. SELEX has 3 lead 
glass calorimeters covering much of the forward solid 
angle. The energy scale for the detector was set 
first by using electron beam scans. Then 7r° decays 
were reconstructed from exclusive trigger data, which 
selected low-multiplicity radiative final states: rj — > 
77 and tt+tt^tt'', oj -^ tt+tt^tt^ as well as 77' and f(1285) 
mesons 16]. The final energy scale corrections were de- 
veloped using 7r° decays from the high-multiplicity charm 
trigger data. Further checks in the charm data set were 
made using single photon decays, e.g., TP -^ A7. The 
uncertainty in the photon e nerg y scale is less than 2%. 
Details can be found in Ref. |l(i| . 

We selected 77 — > 77 candidates in the 77 mass 
range 400-800 MeV/c^. Each photon of the pair has 
E-y >2 GeV. The photon pair has E^^ >15 GeV. The 77 
mass distribution from 10^ charm-trigger events (0.1% of 
the data) is shown in Fig. ^ where a rj signal over a large 
combinatoric background is seen. A fit to a Gaussian plus 
an exponentially falling background yields an rj mass of 
544.8 ± 2.9 MeV/c^ consistent with the PDG value [ll ]. 
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FIG. 1: M(77) distribution for photon pairs in the 77 mass 
region from 0.1% of the data sample. Results for the fit shown 
are in Table Q] The inset shows the background subtracted 77 
signal. The dark points indicate the 77 signal region. 



The mass uncertainty for this and all subsequent states 
is only statistical. 

The observed resolution is 28 ± 4 MeV/c^ consistent 
with the SELEX simulation resuh, 30 MeV/c^. The 
simulation includes all the material in the spectrometer 
and also reproduces the observed tt'^ width as a func- 
tion of energy 46]. The rj signal fraction is 5% within 
a ±60MeV/c^ mass window when we eliminate events 
with more than 5 77 candidates in this mass window. 

We searched for high-mass charm-strange decays that 
followed the pattern Dg plus pseudoscalar meson. We 
had good acceptance and efficiency for the Dg 77 chan- 
nel. The event selection used included the eta selection 
above, and the Dg selection described above, which yields 
a S/N of 4/1. The Ds momenta are typically 150GeV/c 
in the SELEX data set; the E^ > 15 GeV energy cut 
is very loose. We rejected events in which there were 
more than 5 rj candidates in the signal region. This cut 
removed 18 Dg candidates (3.3%) while reducing the 77 
candidate list by 20%. The rj signal region is shown in 
Fig. 2] The final sample consisted of 615 77 candidates 
from 526 Dg candidates. 

The results of our search are shown in the M(KK7r* 77) 
- M(KK7r*) mass difference distribution in Fig. Ufa) . In 
this plot we fixed the 77 mass at the PDG value [13 by 
defining an 77 4-vector with the measured rj momentum 
and the PDG rj mass. A clear peak is seen at a mass 
difference of 666.9 ± 3.3 MeV/c^. 

To estimate the combinatoric background, we matched 
each Dg candidate with 77 candidates from 25 other sam- 
ple events to form a event-mixed sample representing the 
combinatoric background of true single charm production 
and real 77 candidates. The event-mixed mass distribu- 
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TABLE I: Fit results for Figures 1-3 and i3*+(2112) ^ D+7 



Mass 2635.4 + 3.2 
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FIG. 2: M(KK7r± 77) - M(KK7r± ) mass difference distribu- 
tion. Charged conjugates are included. The darker shaded 
region is the event excess used in the estimation of signal sig- 
nificance. The lighter shaded region is the event-mixed com- 
binatoric background described in the text. The inset shows 
the difference of the two with a Gaussian fit to the signal. 
Results for the fit shown are in Tabled 



tion was then scaled down by 1/25 to predict the coni- 
binatoric background in the signal channel. As can be 
seen in Fig. |21 the event-mixed background models the 
background shape very well, but produces no signal peak. 

To estimate the signal yield we subtracted the combi- 
natoric background (light shaded area) from the signal 
data. The resulting difference histogram is plotted in the 
inset in Fig. |5] in the mass-difference range appropriate 
to our search (D^ masses up to 2900 MeV/c^). Out- 
side the peak region the data scatter about 0. The width 
in the D'^K+ mode, to be discussed below, is consistent 
with a 4.9 MeV/c^ Gaussian. We are insensitive to the 
natural width in the D+77 mode. Therefore we fit the dif- 
ference histogram with only a Gaussian with no residual 
background terms. The Gaussian width is fixed at the 
simulation value of 10.9 MeV/c^. The fit yield is 43.4 
± 9.1 events at a mass of 2635.4 ± 3.3 MeV/c^. The 
reduced x^ is 1-10 with a confidence level of 31%. 



To assess the Poisson fluctuation probability to ob- 
serve this excess we note that the resolution is about 1 
bin; we take a 6-bin cluster as a signal region and per- 
form a counting experiment. There are 101 events over 
a background of 54.9 ± 1.5 events giving an excess of 
46.1 events in 6 adjacent bins. The Poisson fluctuation 
probability for this excess is 3 x 10~^ including the un- 
certainty in the background. A conservative estimate of 
the fluctuation probability anywhere in the search region 
(up to 2900 MeV/c2) is 6 x lO"*^. 

The signal does not change with variations of ± 2% in 
the photon energy scale. We also studied combinations 
of events in the Dg mass sidebands with rj candidates and 
candidates in the Dg mass peak with events in the 77 mass 
sidebands. In all cases only smooth combinatoric back- 
grounds, as in Fig. |21 were observed. The rj population 
in this sample, corrected for 77's from the signal channel, 
has a signal fraction of 0.12 ± 0.05, 1.40- higher than the 
global average of 0.05. It is not possible to separate sta- 
tistical fluctuation from a possible 77 enrichment of the 
overall Dg sample. 

A GEANT simulation was also used to determine the 
overall acceptance for these signals. If we detected the 
Ds from a D+ (2632) -^ D+77 decay, then 35 ± 2% of the 
time we also detected the 77 ^ 77. This acceptance was 
obtained by embedding D J (2632) -^ D+77 decay events 
in existing events from the real data. About 55% of 
the Dg decays in SELEX come through this high mass 
state. For comparison, about 25% of the Dg come 
from Dg (2112) decays and a small fraction from either 
DJj(2317) or 0+1(2460) decays, which are marginally vis- 
ible in SELEX data. 

The decay D+(2632) -* D'^K+ is kinematically al- 
lowed. After finding the D^(2632) -^ ^tv signal we 
searched for this second decay mode as confirmation. The 
Z?" sample is the S^ induced Z?" — > if~7r+ subset of the 
sample used in our measurement of the D^ lifetime [13 
with tight D° cuts {L/a > 6, point-back x^ < 5^ ^md 
a good fit to the secondary vertex; x^/i^d < 3). The 
K+ track is > 46 GeV/c (RICH kaon threshold) and is 
strongly identifled by the RICH as > 10 times more likely 
to be a kaon than any other hypothesis. 

The results are shown in Fig. Of a) where we see both 
the known D+j(2573) state clearly and another peak 



above the D^(2573). We fit each peak with a Breit- 
Wigner convolved with a fixed width Gaussian plus a 
constant background term (as suggested from the wrong- 
sign data discussed below). The Gaussian resolution is 
set to the simulation value of 4.9 MeV/c^. The mass dif- 
ference and width of the D^(2573) returned by the fit, 
AM = 705.4± 4.3 MeV/c^ and T = U^l MeV/c^ , re- 
spectively, agree well with the PDG values [Tg] of AM = 
707.9 ± 1.5 MeV/c2 and T = ISt^MeV/c^. The fitted 
mass difference of the second Breit-Wigner is 767.0 ± 2.0 
MeV/c^, leading to a mass for the new peak of 2631.5 ± 
2.0 MeV/c2. The fitted yield is 13.2 ± 4.9 events. The 
signal spread is consistent with the Gaussian resolution, 
even when plotted in 2.5 MeV/c^ bins, limiting the pos- 
sible natural width. For the Breit-Wigner fit we find 
a hmit for the width of < 17 MeV/c^ at 90% confidence 
level. This signal has a significance {[S—B]/^B) of 5.3 a 
in a ±15 MeV/c^ interval. The mass difference between 
this signal and the one seen in the D+77 mode is 3.9 ± 
3.8 MeV/c^ , statistically consistent with being the same 
mass. Unlike the Ds case, the D^K"*" decay contributes a 
small fraction to the SELEX D" sample. 

Combinatoric background will be equally likely to pro- 
duce a D°K~ combination (wrong-sign kaon) as a D"K+. 
The wrong-sign combinations are shown in Fig. |2Ib). 
There is no structure in these data, which fits well to 
a constant background. We conclude that the peak at 
2631.5 MeV/c^ is real and confirms the observation in 
the D+7y mode. 

The relative branching ratio r(D*'K+)/r(D+77 ) must 
be corrected for the relative acceptances of D°, Dg, 
rj, and K~^ mesons, for the r], D^ and Df branching 
ratios, the relative acceptances of the D^(2632) final 
states. We estimate the relative acceptance ratio from 
simulation as 91 ± 3%. The relative branching ratio 
is r(D"K+)/r(D+r/ ) = 0.14 ± 0.06. Relative phase 
space favors the D°K+mode by a factor of 1.53, making 
this low branching ratio even more surprising. 

In conclusion we combined our clean sample of 
Dg mesons with additional photon pairs which made rj 
candidates to study the Dgij mass spectrum. We observe 
a clear peak of 43.4 ±9.1 events with a Poisson fluctua- 
tion probability < 6 x 10~^ at a mass difference of 666.9 
± 3.3 MeV/c^ above the ground state Dg. The back- 
ground shape is well-represented by combinatoric back- 
ground from event-mixing, as discussed above. A cor- 
responding mass peak is also seen in the D^K+channel 
with a significance of 5.3 a at the same mass. The com- 
bined measurement of the mass of this state is 2632.5 ± 
1.7 MeV/c^. The SELEX mass scale systematic error is 
under study, but all charm meson masses measured in 
SELEX agree with PDG values within 1 MeV/c^. The 
state is very narrow, consistent with a width due only to 
resolution in the D°K+decay mode. The 90% confidence 
level upper limit for the width is < 17 MeV/c^. 

SELEX reports these peaks as the first observation of 
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FIG. 3: (a) Ds (2632) -* D°K+ mass difference distribution. 
Charged conjugates are included. The shaded regions are the 
event excesses used in the estimation of signal significances. 
Results for the fit shown are in Tabic |I] (b) Wrong sign back- 
ground D'' K~ events, as described in the text. 



yet another narrow, high-mass Dg state decaying strongly 
to a ground state charm plus a pseudoscalar meson. The 
mechanism which keeps this state narrow is unclear. The 
D"K+ channel is well above threshold, with a Q value 
^ 275 MeV. The branching ratios for this state are also 
unusual. The D+?7 decay rate dominates the D^K+rate 
by a factor of '--^ 7 despite having half the phase space. 
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Note Added in Proof - At the summer conferences the 
B-factory experime nts [lil 1201 and the FOCUS photopro- 
duction experiment [2l| all reported negative results in 
their search for Ds.]{2Q'i2) production. The SELEX data 
have an unusually low ratio of Rs = r'*"'"(2112)/£'+ pro- 
duction (Table I) compared to e+e^ machines. The SE- 
LEX Rs ratio can be brought into agreement with the 
published CLEO result |2| if we exclude the Dg events 
that come through D^j (2632) decays (55% of the SELEX 
Ds yield). This fact, together with the strong production 
asymmetry in the SELEX data, can be interpreted as 
two-component production of charm-strange states from 
the E~ beam. One component involves normal fragmen- 
tation like in photoproduction, the other involves a dif- 
ferent mechanism connected with the beam hadron. To 
understand this production conundrum and to place this 
new state in the spectroscopy of the charm-strange me- 
son system will require careful study from a number of 
experiments in the future. 
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